Well-defined quasiparticles in interacting metallic grains 
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We analyze spectral functions of mesoscopic systems with large dimensionless conductance, which 
can be described by a universal Hamiltonian. We show that an important class of spectral functions 
are dominated by one single state only, which implies the existence of well-defined (i.e. infinite- 
lifetime) quasiparticles. Furthermore, the dominance of a single state enables us to calculate zero- 
temperature spectral functions with high accuracy using the density-matrix renormalization group. 
We illustrate the use of this method by calculating the tunneling density of states of metallic grains 
and the magnetic response of mesoscopic rings. 



The pairing Hamiltonian of Bardeen, Cooper and 
Schrieffer (BCS) is established as the paradigmatic 
framework for describing superconductivity 0, U The 
BCS solution is, however, an approximate one, valid (and 
exceedingly successful) only as long as the mean level 
spacing d is much smaller than the superconducting band 
gap Abcs E| • One of the main features of the BCS 
solution is the description of the excitation spectrum by 
well-defined (i.e. infinite-lifetime) Bogoliubov quasiparti- 
cles, responsible for many of the features of the super- 
conducting state. 

In this Letter, we address the question whether this 
quasiparticle picture prevails in the entire regime of pa- 
rameters - including the case that the samples are so 
small or so weakly interacting that the BCS solution is 
inapplicable - by analyzing spectral functions. For exam- 
ple, the spectral function corresponding to the (noninter- 
acting) particle creation operator ct is given, within the 
BCS solution, by a sharp line in fc-w-space; this reflects 
the infinite lifetime of the quasiparticles. For an inter- 
acting system, this is a very peculiar property, since the 
interactions usually shift a significant portion of the spec- 
tral weight to a background of excitations, responsible 
for the finite lifetime of the quasiparticles. Here we show 
that the unusual property of finding well-defined quasi- 
particles persists to a very good approximation over the 
entire parameter range of the pairing Hamiltonian, and 
is not merely a property of the mean field approximation 
in the BCS regime. We also give a condition for more 
general spectral functions to show analogous behaviour. 

Of central importance is that this result is relevant not 
only in the context of mesoscopic superconductivity, but 
more generally for disordered systems with large dimen- 
sionless conductance g (defined as the ratio between the 
Thouless energy and the mean level spacing d). This is 
because to lowest order in g , the electron-electron in- 
teractions can be described by a remarkably simple uni- 
versal Hamiltonian (UH) which has, besides the 
kinetic energy term H° = ^2, ia e i c lcr c icr) on ly three cou- 



plings: 

H c = E c h 2 , H S = J S S 2 , H p - 



-Xd 



(1) 



Here, E c , J s and —Xd are coupling constants. The sum 
includes all energy levels up to some cutoff lo co at the 
Thouless energy, denoted by the set M. It turns out that 
H c and H s do not affect our result, because they com- 
mute with H° + H p and thus leave the eigenstates invari- 
ant. Therefore, it suffices to take H p - the BCS pairing 
Hamiltonian - as the only interaction term. Therefore, 
for our purposes the difference between the BCS model 
and the UH is only in the cutoff uj co , being at the Debye 
energy for the former and at the Thouless energy for the 
latter. In any case, we define Abcs = w co / sinh(l/A). 

The fact that the zero-temperature spectral function 
Aq (ui) of an operator O is sharply peaked translates to a 
strong condition on the matrix elements of the Lehmann 
representation, which is given by 

A 6 {u) = Y,{&\&W\6\&)S& ~ Ej). (2) 

1-0 

Here |gs) denotes the ground state, 1 1) the excited states 
with energies Ej . For only one sharp peak to be present 
in the spectral function, the sum in Eq. (J2J) must be dom- 
inated by one single eigenstate, say |/)°, whereas all other 
states | J) 7^ |J)° do not contribute. Obviously, it will de- 
pend on the operator O whether this is the case, and if 
so, which is the state |7)°. We show that it suffices that 
O satisfies a rather unrestrictive condition, given after 
Eq. © below and fulfilled for many physically relevant 
quantities. Furthermore, we show that under this con- 
dition, the state |7)° is from a very limited subset of all 
possible excitations, which we characterize below as the 
"No-Gaudino states" . Our finding of well-defined quasi- 
particles therefore implies that only these No-Gaudino 
states are relevant for many physical properties of sys- 
tems that satisfy the conditions of the UH. 

Calculating the spectral function, Eq. (J2J), is usually a 
formidable task, equivalent to diagonalizing the Hamil- 
tonian. Although an exact solution 0, Q exists for the 
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Hamiltonian H p , its complexity in practice does not al- 
low to calculate spectral functions from it. Instead, we 
use the density-matrix renormalization group (DMRG) 
method 8] for this purpose, a numerical variational ap- 
proach that has already be en p roven very useful for an- 
alyzing this model jjj Hoi For suitable operators 
O, we are able to obtain the spectral function from the 
DMRG without the usual complications 0, 0] , because 
the state |J}° - the only one that contributes significantly 
to the spectral function - can be constructed explicitly. 
The existence of a sum rule allows us to quantify the 
contribution of other states \I) ^ \lf, which we find to 
be negligibly small. Finally, we illustrate the use of our 
method of calculating spectral functions by evaluating 
the tunneling density of states and the magnetic response 
of mesoscopic rings. 

Excitation spectrum and No-Gaudino states: Let us 
begin by describing the excitations of the Hamiltonian 
H' p in Eq. 0J. H p has the well-known property that 
singly occupied energy levels do not participate in pair 
scattering; hence their labels (and spins) are good quan- 
tum numbers. Therefore, all eigenstates for which some 
levels i are singly occupied are (as far as the remain- 
ing levels are concerned) identical to those of a system 
with TV in Eq. replaced by Af\B, where B is the set of 
singly occupied levels i |j| . A given state can thus contain 
two kinds of excitations: Pair-breaking excitations that 
go hand in hand with a change of the quantum numbers 
yB, and other many-body excitations that do not. The 
latter were studied in [14j and dubbed "Gaudinos". In 
this spirit, we define the No-Gaudino state as the lowest- 
energy state within a certain sector of the Hilbert space 
characterized by the quantum numbers B. As is shown 
below, this state is easily obtained within the DMRG 
algorithm. 

Let us now specify under which condition the spectral 
function, Eq. (j2J, is dominated by such a No-Gaudino 
state. Any operator can be written as a linear superpo- 
sition of operators 



O 



■ c ■ c 



(3) 



Creating linear superpositions poses no difficulties what- 
soever, therefore it is sufficient to consider operators of 
this form. The central condition we impose on O is that 
all indices ix, ■ • • ,ji be mutually different. O then takes a 
state with no singly occupied levels, B = {}, to the sector 
of the Hilbert space characterized by£> = {ii, ■ ■ ■ , We 
show below that under the above condition, O moreover 
has the crucial property that when acting on the ground 
state, it creates to an excellent approximation the No- 
Gaudino state in this sector. Therefore, the state 0|gs) 
contributing to the spectral function, Eq. J5J), is seen to 
be not only a well-defined eigenstate of the system, but 
moreover a No-Gaudino state. 

In the BCS limit d<A B cs (i-e. at A>l/lnA, where 
N is the number of energy levels within w co ), this follows 



from the identity 

0\gs) =v ix •■•v ik u h ■■■u n \i~' Jl (4) 
^■■■h l ) Q =lt { ^y--lU\9s), (5) 

where the state in Eq. (JjjJ is the No-Gaudino state. Here, 
u, v and 7 are the coherence factors and the Bogoli- 
ubov quasiparticlc operators from BCS theory, as defined 
e.g. in 0. 

In the opposite limit Abcs ^ d (A <C 1/ In N) , where 
perturbation theory in A is valid [ij, the same conclusion 
is obtained: to first order (i.e. up to errors of order A 2 ), 
(5|gs) again creates precisely the No-Gaudino state. 

There is no such simple analytic argument that the 
Gaudino admixture to 0|gs) in Eq. J2J will be negligible 
also in the intermediate regime. However, this assertion 
can be checked numerically by a sum rule, which follows 
from Eq. J5J : 



A(u)<b> = J2(gs\&\I)(I\d\gs) = (gs\&6\gs). (6) 



We define the lost spectral weight wl = (gs\0^0\gs) — 
\(gs\O^\I) \ 2 as the part of Eq. © that is not carried 
by the No-Gaudino state |7)°, but instead lost to other 
background states. As is shown in Fig.^below, this lost 
weight turns out to be negligibly small. 

DMRG algorithm: We now give a brief description of 
the DMRG algorithm as applie d to the UH; more details 
are described elsewhere |l0j, . Energy levels are added 
one by one to the system until it obtains its final size. 
For simplicity, we assume the energy levels to be equally 
spaced, although none of our methods require this as- 
sumption. After adding a level, only a limited number 
m of basis vectors are kept, such that the size of the 
Hilbert space remains numerically manageable. These 
basis vectors are selected in order to represent a num- 
ber of so-called target states accurately; this is achieved 
by the DMRG projection described in By varying 
m between 60 and 140, we estimate the relative error in 
the spectral function from the DMRG projection to be 
of the order of ~ 10~ 5 (for m = 60). This accuracy can 
be improved by increasing m. 

In order to calculate the spectral function correspond- 
ing to the operator O in Eq. ®, we use as target states 
the ground state and a state representing the No-Gaudino 
state \ii ai v ■ ■ ;jf')°, in the BCS limit given by Eq. J5J. 
The sum rule, i.e. the rhs. of Eq. ©, is evaluated in a 
separate run with |gs) and O^gs) as the target states. 

Dominance of a single No-Gaudino state: The fact 
that the spectral function is dominated by one single No- 
Gaudino state is displayed in Fig. ^ Here, the expecta- 
tion value (gs I c i(T cj^ | gs) , which occurs in the sum rule, 
Eq. ijfjjl, with O — q, is plotted (for i — 10, i.e. 10 levels 
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FIG. 1: The matrix element {ci a c\ a ) from Eq. © with O = 
c| CT (dashed line) and the contribution from the No-Gaudino 
state (solid line) as a function of A. Here, i — 10 levels above 
i?Permi out of a total of N — 2 • 40 energy levels. The lost 
weight iul, i.e. the difference between both, is plotted in the 
inset. It shows a maximum in the intermediate regime where 
d ~ Abcs, but even there, iul is less than 0.2% of the total 
spectral weight. 



FIG. 2: The tunneling density of states v(E) for Abcs = 
0.07cj co and N = 10,60,120,400 energy levels in the No- 
Gaudino approximation (solid line; for the sake of better visi- 
bility, the delta peaks in Eq. JS) have been replaced by Gaus- 
sians of width 0.005oj co ). The familiar gap Abcs emerges 
during the crossover from the few-electron (d S> Abcs) to 
the bulk limit (d <C Abcs). In the latter limit, we observe 
agreement with the BCS result (dashed line). 



above £^ e rnii) against the coupling A. It is practically in- 
distinguishable from the contribution | °{^| cJ CT | gs) | 2 from 
the No-Gaudino state only. 

The lost weight wl, shown in the inset of Fig. ^ is 
seen to be less than 0.2% of the total spectral weight 
throughout the entire parameter regime (for i = 10; the 
plots for other values of i, not shown, look similar. The 
maximum lost weight somewhat increases as the level i 
approaches i?Fcrmi, but always remains below 1% of the 
total weight). The lost weight is seen to be vanishingly 
small for small A, as expected in the perturbative regime 
A<Cl/ln7V. Interestingly, the lost weight also decreases 
for large A. This is very untypical for interacting sys- 
tems, and the underlying reason is that the dominance 
of the No-Gaudino state is protected also in the BCS 
regime A> 1/lniV, see Eq. Q. Consequently, the lost 
weight displays a maximum in the crossover regime at 
A ~ 1/lniV. Not shown: We confirmed numerically that 
the coupling A max (^V), at which the lost weight reaches 
its maximum, always scales linearly with 1/ In N as ex- 
pected. The maximum value uiL(A max (iV), N) turns out 
to be a monotonically decreasing function of N. 

Applications: The dominance of the No-Gaudino state 
in the spectral function is not only remarkable by itself, 
but has also high practical value: it allows us to calcu- 
late the spectral function with high precision using the 
DMRG in what we call the "No-Gaudino approximation" 
(NGA), in which only the No-Gaudino state is kept in 
Eq. (J2J). From the spectral function, in turn, many im- 
portant physical quantities can be obtained. The lost 
weight wl, defined after Eq. ©, controls the quality of 
this approximation: when wl vanishes, the NGA is exact. 



As a simple first application, we calculate the tunnel- 
ing density of states v{uo) — 5^ fo ..A.t ( w ) (for w > 0). 
Fig. |21 illustrates that during the crossover from the few- 
electron (d 3> Abcs) to the bulk limit (d <C Abcs), the 
familiar BCS gap of width Abcs emerges together with 
a strongly pronounced peak at oj ~ Abcs as the quasi- 
particle energies are kept away from the Fermi surface 
by the pairing interaction and accumulate at Abcs • The 
lost weight is found to never exceed fractions of 1%, thus 
confirming the accuracy of the NGA. 

As a second example of a quantity that is well cap- 
tured by the NGA, we calculate the prediction of the 
pairing Hamiltonian for the magnetic response AE 2 of 
small metallic rings |l5|. i.e. the derivative of the per- 
sistent current with respect to flux at zero flux. For 
A<Cl/lniV, AE 2 was calculated in perturbation theory, 
for A 3> 1/lnJV the BCS approximation was used |18j . 
Here, using the NGA, we calculate AE 2 for all values of 
A > 0, and specifically in the crossover regime between 
the perturbative and the BCS regimes. 

The linear response to the magnetic flux through a ring 
is given by AE 2 = E% aT 



E 2 ia , where 



£f ar = - 



mn,\I) 



\Pn 



\(I\cl 



ml nf 



il C ml\9s)\ 2 



El 



ar (7) 
and i?2 ia equals the A = value of £"| ar |l8|. Here, m e ,e 
are the electron mass and charge, L is the circumference 
of the ring. P mn is the momentum operator between the 
disordered 1-particle states, labelled by m and n. In the 
highly diffusive regime (to co < 1/r, where r is the elastic 
mean free time), which we assume here for simplicity, 
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FIG. 3: Upper part: The magnetic response AS2 from 
Eq. J7J as function of A for N = 80 energy levels. The upper 
and lower bound from the No-Gaudino approximation (solid 
and dashed line) practically coincide and interpolate between 
the perturbative and the BCS results (dashed-dotted and dot- 
ted line), valid for small and large A, respectively. 
Lower part: The ratio between the No-Gaudino approxima- 
tion and the other results. In the regime A ~ 1/lnY, the 
perturbative and the BCS results underestimate the true re- 
sult (which must lie between the upper and lower bound) by 
a factor of more than 2. 



\P m n\ 2 = P 2 can be taken to be constant for m ^ n, and 
zero otherwise . E\ ar can then be extracted from the 
spectral functions for O mna = c] na c n(J . 

In the NGA, only the states \i — m,n;j — n,m)° 
are retained in Eq. (JJJ . Because the contribution of the 
other states, which are neglected, is always positive, the 
NGA produces a lower bound for Ei. An upper bound 
can be found as well, namely by replacing the energy 
denominator of Eq. Q by the energy £^ n of the No- 
Gaudino state, which is known to be smaller than the 
energy of any other contributing state. Then, the sum 
over 1 7) can be eliminated, and the resulting expression 
for the upper bound is 



E>=a 



E 



(gs\(cl 



-C' 1 C 1 )fcl 



"mi 1 



' C ni C mi 



)\gs) 



E° 



(8) 

The results of our calculation are presented in Fig. |3 
where the upper and lower bound is compared to the per- 
turbative and to the BCS result, given in [lq. The lower 
and upper bounds practically coincide (with an error of 
< 0.5%) in the entire parameter regime; this reflects the 
high accuracy of the NGA. As expected, the perturbative 
result is reproduced for small A (Abcs ^ d), the BCS 
result for large A (Abcs ~> d). However, both results un- 
derestimate the exact result by a factor of up to 2.5 in a 
large intermediate regime (Fig. [3] bottom). Interestingly 
we find that the magnetic response is much larger than 



the BCS value also in a regime in which A ^> d, where the 
BCS approximation is expected to be valid. This is due 
to a large contribution of the distant levels from A B cs U P 
to the interaction cutoff, which the BCS approximation 
neglects. A similarly large contribution from distant lev- 
els has previously been found also for the condensation 
energy 0] and single particle properties [l9|. 
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